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Abstract 

A  voltage  degradation  model  for  the  low-pressure  proton  exchange  membrane  fuel  cell  (PEMFC)  stack  used  in  a  fuel  cell  bus  is  presented: 
(1)  the  oxygen  concentration  term  was  derived  from  the  PEMFC  output  voltage  equation,  and  the  concept  of  oxygen  concentration  resistance 
coefficient  was  introduced;  (2)  a  5  kW  low-pressure  PEMFC  stack  was  used  in  this  study.  Two  similar  tests  were  carried  out  before  and  after  the 
stack  operating  in  the  driving  cycle  for  640  h.  First,  the  ohmic  losses  under  different  temperatures  were  measured  using  the  current  interrupt  method 
and  formulized  with  linear  fitting  method.  Then,  the  oxygen  concentration  term  was  studied  by  the  experiments  with  different  air  stoichiometric 
ratios  while  keeping  the  other  operating  parameters  unchanged.  The  oxygen  concentration  resistance  coefficient  was  obtained  from  the  difference 
of  voltages  for  the  PEMFC  stack  in  different  air  stoichiometric  ratios  using  the  genetic  optimization  algorithm.  Then,  the  activation  loss  was 
obtained  based  on  the  PEMFC  output  voltage,  the  ohmic  loss,  and  the  concentration  loss.  The  degradation  model  of  the  stack  was  built  finally  by 
comparing  the  two  test  results;  (3)  the  correlation  of  the  model  to  the  actual  experimental  data  is  good;  (4)  the  overvoltage  of  the  stack  with  aging 
was  analyzed  using  this  model.  The  analysis  showed  that  the  activation  overvoltage  dominated  the  stack  loss  with  about  80%  of  the  total  losses, 
followed  by  the  ohmic  loss.  The  concentration  loss  almost  does  not  change  with  aging  in  the  driving  cycle  condition;  (5)  the  comparison  of  the 
simulation  with  the  actual  data  from  the  PEMFC  bus  running  for  30,000  km  indicated  that  after  36,000  km  the  rated  power  of  the  PEMFC  bus  must 
be  reduced. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

As  the  vehicle  emission  regulations  become  more  and  more 
strict,  electrical  vehicles  have  become  the  focus  of  people’s 
attention.  Proton  exchange  membrane  fuel  cell  (PEMFC)  is 
one  of  the  most  promising  power  sources  for  electric  vehicles. 
PEMFC  stack  is  the  core  of  the  PEMFC  vehicle.  The  perfor¬ 
mance  and  the  lifetime  of  the  fuel  cell  system  and  the  fuel 
cell  vehicle  are  strongly  influenced  by  the  degradation  of  the 
PEMFC  stack.  Therefore,  it  is  necessary  to  establish  a  voltage 
degradation  model  of  the  PEMFC  stack.  There  are  two  type  of 
models  [1],  one  is  the  mechanism  model,  and  the  other  is  the 
empirical  model.  The  mechanism  model  is  established  based  on 
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some  proper  assumptions.  It  consists  of  the  mass  transfer  func¬ 
tions,  the  thermal  transport  functions,  and  the  electrochemical 
reaction  functions.  Mechanism  model  is  complex  though  it  can 
discover  the  inner  mechanism  of  the  fuel  cell  stack.  It  is  more 
suitable  for  studying  the  inner  phenomena  of  fuel  cell  stacks. 
The  empirical  model  described  with  meaningful  coefficients  is 
usually  based  on  the  fuel  cell  voltage-ampere  experimental  data 
without  considering  the  fuel  cell  inner  structure.  The  empirical 
model  is  much  simpler  and  more  suitable  for  the  system  analy¬ 
sis.  Different  empirical  models  have  been  developed  [2-14].  The 
model  developed  by  Amphlett  et  al.  [6-9]  uses  a  combination 
of  mechanistic  and  empirical  modeling  techniques,  and  every 
coefficient  has  its  meaning.  Based  on  the  Amphlett’ s  model, 
Mann  et  al.  took  into  the  reactant  humidity  in  their  study  [1 1].  A 
more  sophisticated  model  has  been  developed  by  Pukerushpan 
et  al.  which  incorporates  the  reaction  gases  humidity  and  the 
concentration  overvoltage  into  the  Amphlett’s  model  [12,13]. 
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Nomenclature 

co2,  interface  oxygen  concentration  at  the  cathode  mem¬ 
brane/gas  interface  (mol  cm-3) 

Cr  oxygen  concentration  resistance  coefficient 

(atm  mA  1  cm2) 

E  thermodynamic  potential  (V) 

i  current  density,  (mA  cm-2) 

in  internal  and  fuel  crossover  equivalent  current  den¬ 

sity  (mAcm-2) 

I  stack  output  current  (A) 

n  cell  number  in  series 

/?anode, channel  the  total  pressure  in  the  anode  gas  channel 
(atm) 

pC2L  the  cathode  operating  pressure  (atm) 

Ph2  , interface  hydrogen  partial  pressure  in  the  reactive  inter¬ 
face  (atm) 

pn2o, anode  the  partial  pressure  of  water  vapor  in  the  flow 
channels 

Po2, interface  oxygen  partial  pressure  in  the  reactive  inter¬ 
face  (atm) 

Po2 , anode  oxygen  partial  pressure  in  the  bipolar  channel 
(atm) 

Po2, channel, a.i  ,  Po2, channel, z2  the  oxygen  partial  pressure  in 
the  bipolar  channel  under  stoichiometric  ratio 
and  A. 2  (atm) 

Panode  h2o  water  saturated  pressure  in  the  anode  (atm) 
PcslU2o  the  water  saturation  partial  pressure  in  cathode 
(atm) 

R  stack  internal  resistance  (£2) 

t  stack  operating  time  under  the  driving  cycle  test 

(h) 

T  fuel  cell  operating  temperature  (K) 

V  fuel  cell  output  voltage  (V) 

Vact  activation  overvoltage  (V) 

Vcon  concentration  overvoltage  (V) 

Vohm  ohmic  overvoltage  (V) 

Vstack  stack  output  voltage  (V) 

A  V  the  voltage  difference  of  the  PEMFC  output  volt¬ 

age  under  the  air  stoichiometric  ratio  X\  and  X2 
(V) 

Greek  letters 

r ]  fuel  cell  stack  efficiency 

§i>§2>§3>£4>fi>f2>f4  empirical  coefficients  for  calcula¬ 
tion  of  activation  overvoltage 
£3,  \jf  empirical  coefficients  for  calculation  of  concen¬ 
tration  overvoltage 


The  PEMFC  performance  becomes  worse  as  the  stack  ages.  All 
the  above  models  deal  little  with  degradation.  A  few  empiri¬ 
cal  models  can  handle  the  voltage  degradation  but  were  used 
for  the  system  analysis,  i.e.,  Fowler  et  al.  have  incorporated  the 
voltage  degradation  into  their  model  based  on  the  Amphlett’s 
model  [14].  In  this  paper,  we  report  a  PEMFC  stack  model  used 


for  both  the  fuel  cell  bus  hybrid  system  and  the  fuel  cell  sys¬ 
tem  analysis.  We  introduce  this  Amphlett’s  model-based  voltage 
degradation  model  for  low-pressure  PEMFC  stack  into  our  fuel 
cell  bus.  We  first  derive  the  oxygen  concentration  term  from 
the  PEMFC  output  voltage  equation  and  introduce  the  concept 
of  oxygen  concentration  resistance  coefficient,  then  use  a  5  kW 
low-pressure  PEMFC  stack  module  for  the  study.  Two  sets  of 
the  similar  tests  are  carried  out  before  the  driving  cycle  test  and 
after  the  stack  operating  in  the  driving  cycle  for  640  h.  Then 
the  model  coefficients  changing  with  time  are  obtained  from  the 
experimental  data.  Finally,  the  degradation  model  is  used  to  ana¬ 
lyze  the  losses  of  the  fuel  cell  stack.  By  comparing  the  model 
simulation  results  with  the  test  data  from  the  fuel  cell  bus  run¬ 
ning  for  30,000  km,  the  degradation  of  the  stack  used  in  our  fuel 
cell  bus  is  evaluated. 


2.  Development  of  PEMFC  model 


The  output  voltage  of  the  fuel  cell  is  defined  as  the  function 
of  the  thermodynamic  potential,  the  activation  overvoltage,  and 
the  ohmic  overvoltage,  as  shown  by  the  following  Eq.  (1) 

V  =  E  +  Vact  +  Vohm  (1) 


where  V,  Vact,  V0hm  and  E  are  the  fuel  cell  output  voltage,  the 
activation  overvoltage,  the  ohmic  overvoltage,  and  the  thermo¬ 
dynamic  potential,  respectively.  E  can  be  described  via  Nernst 
equation  as 


E  =  1.229  -  0.85  x  10_3(T  -  298.15)  +  4.309  x  10"57 


ln(PH2, interface)  +  ^  ^(/?02,  interface) 


(2) 


Eq.  (2)  can  be  reformed  as 

E  =  E'  +  2.15425  X  10~5r  ln(/?o2, interface)  O) 


E'  =  1.229  -  0.85  x  10_3(r  -  298.15)  +  4.309 

X  10  T  ln(/>H2, interface)  13 ft) 

where  T  is  the  fuel  cell  operating  temperature.  po2, interface  is 
the  oxygen  partial  pressure  in  the  reactive  interface.  /?h2, interface , 
the  hydrogen  partial  pressure  in  the  reactive  interface  can  be 
expressed  as 

PH2,  interface  =  P  anode,  channel  —  PH20, anode  (3b) 

where  Panode, channel  is  the  total  pressure  in  the  anode  gas  chan¬ 
nel.  The  partial  pressure  of  water  vapor  in  the  flow  channels, 
Ph2o, anode  can  be  defined  as  following  Eq.  [6] 

PH20, anode  =  0-5/4node  ,H20  (3c) 

where  ps^odQ  H2q  is  the  water  saturated  pressure  in  the  anode  at 
the  fuel  cell  operating  temperature. 

Activation  overvoltage  can  be  given  as  the  following  Eq.  [6] 

Vact  =  H”  &  T  H”  §3 ^[ln(co2, interface)]  H”  ?4T[ln(/  +  /n)]  (4) 
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where  co2, interface  is  the  concentration  of  dissolved  oxygen  at 
the  gas/liquid  interface,  which  can  be  defined  by  Henry’s  law 
expression  of  the  form  [6] 


c02,  interface 


PO2,  interface 

5.08  x  106exp(— 498/r) 


(4a) 


Supposing  each  cell  performance  in  the  stack  is  the  same,  the 
output  voltage  of  the  stack  can  be  expressed  as 

Vstack=^V  =  1t[E  +  Vact  +  Vohm  +  tyT  ln(/?o2,  channel  —  Gr/)] 

(ID 


§1,  $2,  $ 3 ,  $4  are  constants.  i  +  in  is  the  total  current  density. 
i  is  the  output  current  density,  and  in  is  the  internal  and  fuel 
crossover  equivalent  current  density  [15]. 

Substituting  Eq.  (4a)  into  Eq.  (4),  it  can  obtain: 

Vact  =  fl  +  &T  +  £3  T  [ln(/?o2,  interface)]  +  &TUn(i  +  /n)] 

Vact  =  Vact  +  £3  T  [In  (/?o2,  interface)] 
where  fi,  £2,  £3,  £4  are  constants. 

Vact  =  Cl  +  KiT  +  ^4^[ln(/  +  in)]  (5  a) 

Ohmic  overvoltage  can  be  described  as 

Vohm  —  —Ri  (6) 

where  R  is  the  internal  electric  resistance. 

Combining  Eqs.  (6),  (5),  (3)  and  (1)  yields 

V  =  E'  +  v^t  +  Vohm  +  2.154  X  10-5rin(po2, interface) 

H”  £3  T  In (po2, interface)  (7) 

This  equation  can  be  simplified  to: 

V  =  E  +  Vact  +  V0hm  In (PO2, interface)  (8) 


where  n  is  the  cell  number  in  series  in  the  stack.  The  fourth  term 
in  the  right  side  of  Eq.  (1 1)  is  defined  as  oxygen  concentration 
different  term.  In  Eq.  (11),  E\  Vact,  and  V0hm  do  not  relate  to 
oxygen  partial  pressure.  In  the  following  section,  tests  are  elab¬ 
orately  designed  to  obtain  each  of  the  constants  in  the  above 
equations. 

3.  Experimental 

Tests  are  carried  out  on  the  test  station  FCATS  G700  man¬ 
ufactured  by  GREENLIGHT  power  technologies  company.  Air 
is  from  the  air  station,  and  hydrogen  is  from  the  compressed 
hydrogen  cylinders. 

A  5  kW  low-pressure  PEMFC  stack  module  is  used  for  the 
study.  The  same  kind  of  PEMFC  is  also  used  in  our  fuel  cell 
buses  (one  of  our  fuel  cell  bus  PEMFC  stack  is  about  72  kW). 
The  5  kW  stack  consists  of  80  cells  in  series.  The  active  area  of 
the  PEMFC  is  280  cm2,  and  the  membrane  is  Nation®  1 12.  The 
basic  parameters  of  the  stack  are  listed  in  Table  1 . 

3.1.  Tests  before  the  driving  cycle  test 


where 

f  =  2.154  x  10“5  +f3  (8a) 

Assume  the  oxygen  partial  pressure  in  the  bipolar  channel 
as  Po2, channel?  and  assume  the  oxygen  partial  pressure  in  the 
reaction  interface  as 

PO2,  interface  =  PO2,  channel  —  Cri  (9) 

where  Cr  is  defined  as  the  oxygen  concentration  resistance  coef¬ 
ficient,  which  is  influenced  by  the  structure  of  gas  distribution 
layer  (GDL)  and  the  liquid  water  contained  in  the  GDL. 
Substituting  Eq.  (9)  into  Eq.  (8)  gives: 

V  =  E  +  Vact  +  Vohm  H”  la (PO2, channel  —  Gr/)  (10) 


Before  the  driving  cycle  test,  a  set  of  tests  are  designed  to  dis¬ 
tinguish  the  empirical  coefficients  from  the  equations  in  Section 
2. 

3.1.1.  Ohmic  resistances  tests 

The  ohmic  resistances  are  determined  by  using  the  current 
interrupt  method  [15].  The  electronic  load  is  adjusted  until  the 
desired  test  current  is  stable,  then  the  load  current  is  switched  off 
rapidly.  The  changed  voltage  is  recorded  to  the  storage  oscillo¬ 
scope.  The  internal  ohmic  resistance  can  be  obtained  by  dividing 
the  changed  voltage  with  the  current. 

The  test  results  of  the  reference  [16]  show  that  the  internal 
resistances  of  very  thin  membranes  do  not  change  with  the 
current  density  if  the  reactive  gases  are  well  humidified.  As 


Table  1 

The  basic  parameters  of  the  PEMFC  stack 

Cells  number  in  series  Active  area  (cm2)  Oxidant  Fuel  Rated  power  (kW)  Rated  operating  current  (A) 

80  280  Air  Pure  hydrogen  5  90 

Rated  operating  Rated  operating  Cathode  rated  air  Anode  rated  Operating  Operating 

temperature  (°C)  pressure  stoichiometric  ratio  stoichiometric  ratio  temperature  (°C)  pressure  (kPa) 

60  Atmospheric  pressure  2.5  1.5  Room  temperature  to  65  Atmospheric  pressure  to  50 


Operating  relative 
humidity  (%) 

Cathode 

stoichiometric  ratio 

Anode 

stoichiometic  ratio 

Max  allowable  different 
pressure  between  anode 
and  cathode  (kPa) 

Max  allowable 
continuous  operating 
current  (A) 

The  allowable 

minimum  start 
temperature  (°C) 

60-100 

2. 1-2.7 

1.3— 1.5 

30 

168 

5 
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Current  Density  (mA  crrv2) 

Fig.  1.  The  internal  resistances  of  the  stack  in  different  current  density  and  air 
dew  point  temperatures. 


Nafion  112  is  used  in  the  5kW  low-pressure  PEMFC  stack, 
the  hydrogen  is  always  humidified  around  90%,  and  the 
humidity  of  the  air  is  always  around  60-100%.  We  design 
an  experiment  to  confirm  the  relationship  between  the  stack 
internal  resistance  and  the  current  density.  The  stack  operating 
temperature  is  kept  at  60  °C.  The  reactive  gases  pressure  is 
atmosphere  pressure.  The  air  stoichiometric  ratio  is  2.3.  The 
hydrogen  stoichiometry  is  1.3,  and  the  dew  point  temperature 
of  the  hydrogen  is  55  °C.  The  test  results  show  that  if  the  dew 
point  of  the  air  is  kept  between  50  °C  and  60  °C,  the  internal 
resistances  of  the  stack  change  very  little  with  the  current 
density,  especially  under  high  current  density  as  shown  in 
Fig.  1. 

Other  experiment  results  show  that  under  good  humidity  con¬ 
dition  the  stack  internal  resistances  do  not  change  with  the  air 
stoichiometric  ratio  below  2.5.  We  also  carried  out  the  experi¬ 
ment  to  find  the  relationship  of  the  stack  internal  resistances  to 
the  stack  operating  temperatures.  In  the  experiment,  the  operat¬ 
ing  pressure  is  atmosphere  pressure,  the  air  stoichiometric  ratio 
is  2.3,  the  hydrogen  stoichiometric  ratio  is  1.3,  reactants  are  both 
saturated.  The  results  are  shown  in  Fig.  2.  The  changed  voltages 
are  linear  to  the  interrupted  currents,  and  the  slopes  of  these  lines 
are  greater  as  the  stack  operating  temperatures  reduces.  These 
mean  that  the  internal  resistances  of  the  stack  have  little  rela¬ 
tion  with  the  currents  and  become  higher  as  the  stack  operating 
temperature  is  lower.  The  relationship  of  the  internal  resistance 
to  the  stack  operating  temperature  is  a  conic  as  shown  in  Fig.  3. 


Fig.  2.  The  changed  voltages  with  different  interrupted  currents  in  different 
stack  operating  temperatures. 


<S  0.0361 - 1 - 1 - 1 - 

£  300  310  320  330  340 

S  Stack  Operating  Temperature  (K) 


Fig.  3.  The  relationship  of  the  stack  internal  resistances  with  the  stack  operating 
temperatures. 


The  stack  resistance  can  be  expressed  as  Eq.  (12) 

Rn  =  5.61  x  1(T6  x  T2  -  3.7248  x  10“3  x  T  +  0.655  (12) 
So  the  ohmic  overvoltage  can  be  expressed  as 
Vohm  =  (7.013  X  10“8  X  T2  -  4.656  x  10“5 

x  T  +  8.188  x  10“3)  x  /  (13) 


3.1.2.  The  air  concentration  term  tests 

The  concentration  term  can  be  easily  separated  from  the  stack 
output  voltage  by  the  following  method. 

The  air  stoichiometric  ratio  changes  from  k\  to  X2,  while  the 
other  operating  parameters  kept  unchanged.  Then  based  on  the 
Eq.  (10),  we  have  the  following  equations: 

V\  =  E  -\-  Vact  +  Vohm  H”  ln(po2  channel, ~  Cri)  (14) 

V2  =  E  +  Vact  +  Vohm  T  ln(/?02, channel, A,2  —  l)  (15) 

Subtracting  these  two  equations  gives 

AV  —  Vrr{ln(/702, channel, A.2  —  CrO  —  In  (PO2,  channel,  X 1  —  CrO } 

(16) 


where  AV  is  the  voltage  difference  of  the  PEMFC  out¬ 
put  voltage  under  the  air  stoichiometric  ratio  k\  and  X2. 
Po2, channel, A.!  and /?o2, channel, a2  are  the  oxygen  partial  pressure  in 
the  bipolar  channel  under  stoichiometric  ratio  A,i  and  X2,  respec¬ 
tively. 

Under  good  humidity  condition,  the  oxygen  partial  pressure 
in  the  bipolar  channel  can  be  expressed  as  Eq.  (17) 


PO2, channel, A,  —  (Pea  Pca^t^o)  x  0.21  x 


(k  ~  1  )/k 
(1  -0.21 /A.) 


/In 


'1  —  0.21/A. 

_  (k  -  1  )/k  _ 


(17) 


where  k  is  the  air  stoichiometric  ratio.  pscf  H2o  *s  l^e  water  sat_ 
uration  partial  pressure  in  cathode.  pca  is  the  cathode  operating 
pressure. 

The  constants  x/r  and  Cr  can  be  easily  obtained  by  the  follow¬ 
ing  designed  experiments.  The  stack  operating  temperatures 
are  kept  at  30  °C,  40  °C,  50  °C,  60  °C,  respectively,  and  the  air 
stoichiometric  ratios  are  changed  with  k  =  1.9,  2.1,  2.3,  2.5,  and 
2.7,  respectively.  The  results  V-I  curves  under  the  60  °C  are 
shown  in  Fig.  4.  The  stack  output  voltage  differences  ftAV  of 
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When  T=  303. 15  K,  then  f  =  2.0  x  10“4,  Cr  =  2.168  x  10“4. 

So  the  constant  xjf  changes  very  little  with  the  operating  tem¬ 
perature.  Substituting  \fr  into  Eq.  (8a),  it  gives  £3  =  1.78  x  10-4. 
The  thermodynamic  potential  is  not  sensitive  to  oxygen  con¬ 
centration  by  comparing  2.154  x  10-5  with  £3.  The  oxygen 
concentration  resistance  coefficient  Cr  increases  a  little  as  the 
operating  temperature  decreases,  which  is  due  to  more  liquid 
water  in  the  GDL  at  low  operating  temperature.  According  to 
the  test  data,  Cr  can  be  expressed  as: 

Cr  =  -3.2  X  10“8  XT2  +  1.92  X  10-5  x  T  -  2.663  x  1CT4 

(18) 


Fig.  4.  The  stack  V-I  curves  under  different  air  stoichiometric  ratios. 

the  experiments  are  shown  in  Figs.  5  and  6.  Then,  the  constants 
x\r  and  Cr  can  be  derived  by  the  genetic  arithmetic  method  [17]. 
The  genetic  arithmetic  method  is  used  to  process  the  difference 
voltages  of  the  PEMFC  stack.  The  results  are  shown  in 
following: 

When  T=  333.15  K,  then  ^  =  2.0  x  1(T4,  Cr=  1.800  x  1(T4. 

When  T=  323.15  K,  then  ifr  =  2.0  x  1(T4,  Cr  =  2.007  x  1(T4. 

When  T=  313.15  K,  then  \fr  =  2.0  x  1(T4,  Cr  =  2.100  x  1(T4. 
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X=2.1-X=l  .9  — ■—  2.3-X^l  .9 
X=2.5-X=l  .9  X  X=2.7-X=l  .9 


Fig.  5.  The  stack  voltage  differences  under  different  air  stoichiometric  ratios 
(1). 


3.1.3.  Activation  overvoltage  calculation 

As  E\  Vohm,  xj/T\n(po2, channel  -  Cr/)and  the  PEMFC  out¬ 
put  voltages  V  in  Eq.  (10)  have  been  obtained,  the  constants 
in  V^ct  can  be  also  obtained  by  using  the  same  genetic  arith¬ 
metic  method.  The  results  show  that  £i  and  £2  are  little  sensitive 
to  the  operating  temperature,  £1  =0,  £2  =  3.273  x  10-4,  while 
£4  and  in  are  strongly  affected  by  the  operating  temperature. 
For  example,  when  T=  333.15  K,  323.15  K,  313.15  K,  303.15  K, 
the  related  two  parameters  are  £4  =  8.106  x  10-5, 9.748  x  10-5, 
1.145  x  10-4,  1.316  x  10-4,  and  zn  =  0.99,  1.603,2.314,3.324. 
Constants  £4  and  the  internal  and  fuel  crossover  equivalent  cur¬ 
rent  density  in  increase  a  little  as  the  operating  temperature 
decreases  (zn  is  typically  in  the  range  of  l-5mAcm-2  [18]). 
They  can  be  expressed  as  following: 

£4  =  -1.686  X  l<r6  X  T  +  6.427  x  10“4  (19) 

in  =  9.93  x  10“4  x  T2  -  0.709  x  T  +  127  (20) 

Substituting  them  into  the  Eq.  (10)  gets 
C  =  3.2725  X  10“4  X  T  +  (-1.6864  x  10“6  x  T  +  6.427 
x  10“4)  x  rin[/+(9.93  x  10“4  x  T2-0J09T  +  127)] 

(21) 
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— X^2.3-X^2.1  — X^2.5-X^2.1  — A—  h=2.1-h=2. 1 

Fig.  6.  The  stack  voltage  differences  under  different  air  stoichiometric  ratios 
(2). 


3.2.  Tests  after  640  h  driving  cycle  test 

The  testing  driving  cycle  is  shown  in  Fig.  7  [19].  This  driv¬ 
ing  cycle  was  derived  from  the  hybrid  fuel  cell  bus  simulation 
model  based  on  the  Chinese  typical  bus  driving  cycle.  There  are 
two  phases  (phase  I  and  phase  II)  in  the  Fig.  7.  The  driving  cycle 
consists  of  four  phase  I  and  one  phase  II.  The  5  kW  low-pressure 
PEMFC  stack  has  been  tested  under  the  driving  cycle  for  640  h. 
The  stack  output  voltage  history  under  the  current  100  A  is 
shown  in  Fig.  8.  The  degradation  of  the  stack  output  voltage  is 
roughly  linear.  The  degraded  rate  is  about  5.8  mVh-1,  and  the 
average  cell  output  voltage’s  degraded  rate  is  about  72.5  jxV  h_1 , 
which  is  worse  than  those  presented  in  other  literatures  [14,20]. 
The  stack  V-I  curves  before  and  after  640  h  driving  cycle  test 
are  compared  in  Fig.  9.  The  voltage  has  dropped  about  4  V  after 
640  h  driving  cycle  test  under  the  current  100  A,  which  means 
that  the  stack  voltage  has  degraded  about  6.9%. 
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Fig.  7.  The  PEMFC  stack  test  driving  cycle. 


Fig.  8.  The  stack  output  voltage  degradation  curve  under  the  current  100  A. 


In  order  to  compare  the  losses  of  the  stack  before  and  after 
the  driving  cycle  test,  the  similar  tests  listed  in  Section  3.1  are 
carried  out  after  640  h  driving  cycle  test. 

The  results  of  the  internal  resistances  tests  under  the  operating 
temperature  333.15  K  are  shown  in  Fig.  10.  After  640  h  of  driv¬ 
ing  cycle  test,  the  internal  resistance  rose  from  3.6956  x  10~2Q 
to  4.98  x  10-2  Q  by  34.7%.  The  increasing  rate  is  about 
2.006  x  10-5  £2h-1.  The  increasing  rate  for  converted  aver¬ 
age  cell  internal  resistance  is  2.508  x  10-7  £2h-1.  Assuming 


Before  the  driving  cycle  test  — ■ —  After  640  hours  driving  cycle  test 


Fig.  9.  Comparison  of  the  V-I  curves  before  and  after  640  h  driving  cycle  test. 


Fig.  10.  The  changed  voltages  in  different  interrupted  currents  after  640  h  driv¬ 
ing  cycle  test. 


the  increase  of  the  internal  resistance  is  linear,  we  find  the 
ohmic  loss 

yohm  =  (7.013  x  10-8  x  T2  -  4.656  x  10~5  x  T  +  8.188 

x  10“3  +  2.508  x  10“7r)  x  I  (22) 

where  t  is  the  stack  operating  time  under  the  driving  cycle  test. 

The  concentration  term  tests  shown  that  after  640  h  running 
of  driving  cycle  test,  the  coefficients  of  the  air  concentration  term 
changed  little.  And  under  the  operating  temperature  333.15  K, 
i/f  =  2.0  x  10-4,  Cr=  1.810  x  10-4. 

Using  the  similar  method  listed  in  Section  3.1.3  gets  fi  =0, 
£2  =  3.253  x  10-4,  £4  =  9.882  x  10-5,  in  =  4.676  under  the  oper¬ 
ating  temperature  of  333. 15  K.  Compared  to  those  before  driving 
cycle  test,  £1  and  £2  have  changed  very  little,  and  the  inter¬ 
nal  and  fuel  crossover  equivalent  current  density  in  have  risen 
about  four  times,  which  means  that  there  are  may  be  some 
hot  pinholes  in  the  membrane.  Assuming  the  increase  of  inter¬ 
nal  and  fuel  crossover  equivalent  current  density  in  is  linear, 
it  gives 

l4t  =  3.2725  x  10"4  x  T  +  (-1.6864  x  10-6  x  T  +  6.427 
x  10“4  +  2.775  x  10“80  x  T  In [7  +  (9.93  x  10“4 
x  T2  -  Q.1Q9T  +  127  +  0.00576f)] 


(23) 
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Fig.  11.  Comparison  of  the  model  calculation  to  the  experimental  value 
(the  operating  pressure  is  1.0-1. 2  atmosphere  pressure.  The  operating 
temperature  is  303.15-333.15  K.  The  relative  humidity  of  anode/cathode 
is  80-100%/60-100%.  The  stoichiometric  ratio  of  anode/cathode  is 
1.1-1. 5/1. 7-2.9). 


4.  Model  validation 

The  correlation  of  the  model  to  20  actual  experimental  data 
under  different  operating  temperatures  and  air  stoichiometric 
ratios  is  shown  in  Fig.  11.  The  comparison  of  the  model  average 
cell  V-I  to  actual  cell  V-I  curves  is  shown  in  Fig.  12.  The  model 
voltage  degradation  rate  under  100  A  current  is  about  7.15%, 
while  the  actual  voltage  degradation  rate  is  about  6.9%  as  shown 
in  Section  3.2.  The  precision  of  the  model  is  about  5%,  which 
means  the  above  modeling  method  is  feasible.  In  the  next  section 
we  will  use  the  model  to  predict  and  analyze  the  degradation  of 
the  stack. 

5.  Model  simulation  and  analysis 

The  V-I  curves  and  the  output  powers  of  the  stack  changed 
with  time  are  shown  in  Figs.  13  and  14,  respectively.  The  test 
conditions  are  shown  in  following:  the  operating  temperature: 


4—  BOL  -  -■  500  Hours  A  1000  Hours 


1 500  Hours  — —  2000  Hours  — • —  2500  Hours 
— I —  3000  Hours 

Fig.  13.  The  simulation  V-I  curves  of  the  fuel  cell  stack  in  different  age. 

T  =  333.15  K,  H2  and  air  stoichiometric  ratios:  A  =  1.5  and  2.5, 
respectively. 

In  order  to  analyze  the  losses  of  the  stack  conveniently,  the 
concentration  loss  is  defined  as  the  following  equation 

K:on  =  —  [tyT  ln(/?o2, channel)  —  lnCP02, channel  —  Cvif\  (24) 

Substituting  this  Eq.  (24)  into  Eq.  (10)  gives 

V=E"  +  V"t  +  Vohm  +  Vcon  (25) 

where  E"  is  defined  as  normal  thermodynamic  potential 

E"  =  E'+  2.15425  x  lf)-5rin( /^.channel)  (25a) 

V"t  is  activation  loss 

^act  =  ^act  +  ft  T  ln(/?02, channel)  (25b) 

Under  the  rated  current  90  A  defined  by  the  manufacturer,  the 
losses  of  the  stack  with  aging  are  shown  in  Table  2. 

At  the  beginning  of  the  driving  cycle,  the  normal  thermody¬ 
namic  potential  is  94.73  V,  and  the  output  voltage  of  the  stack  is 
58.21  V.  The  ohmic  loss  is  3.31  V.  The  activation  loss  is  30.39  V, 
and  the  concentration  loss  is  2.80  V,  the  efficiency  of  the  stack 
is  77  =  61 .45%. 


Current  Density(mA  •  cnr2) 


Fig.  12.  Comparison  of  the  model  V-I  curves  to  the  experimental  V-I  curves 
(the  operating  pressure  is  atmosphere  pressure.  The  dew  point  of  anode/cathode 
is  330.15  K/325.15  K.  The  stoichiometric  ratio  of  anode/cathode  is  1.2/2. 3). 


Current  Density  (mA  cnr2) 

♦  -BOL  — ■ — 500  Hours  — A — 1000  Hours  1500  Hours 

— X —  2000  Hours  — • —  2500  Hours  - 1 - 3000  Hours 

Fig.  14.  The  simulation  output  powers  of  the  fuel  cell  stack  in  different  age. 
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Table  2 

The  performance  of  the  stack  with  age  under  the  rated  current  90  A 
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Age  (h) 

Output  voltage  (V) 

Stack  efficiency  (%) 

Activation  loss  (V) 

Ohmic  loss  (V) 

Concentration  loss  (V) 

Total  loss  (V) 

BOL 

58.21 

61.45 

30.39 

3.31 

2.80 

36.5 

500 

55.15 

58.22 

32.55 

4.22 

2.80 

39.57 

1000 

52.08 

54.98 

34.71 

5.12 

2.80 

42.63 

1500 

49.00 

51.73 

36.88 

6.02 

2.80 

45.7 

2000 

45.93 

48.48 

39.06 

6.93 

2.80 

48.79 

2500 

42.84 

45.22 

41.24 

7.83 

2.80 

51.87 

After  2500  h  driving  cycle-test,  the  output  voltage  of  the  stack 
becomes  42.84  V.  The  ohmic  loss  is  7.83  V.  The  activation  loss  is 
41 .24  V.  The  concentration  loss  is  still  2.80  V,  and  the  efficiency 
of  the  stack  is  f)  =  45.22%. 

The  activation  loss  dominates  the  total  losses  in  the  stack 
with  almost  80%.  The  ohmic  loss  is  the  second  largest  loss. 
The  concentration  loss  is  the  mildest  and  does  not  change  with 
aging  in  the  driving  cycle  condition.  After  2500  h  driving  cycle 
test,  the  total  losses  change  from  36.5  V  to  51.87  V  with  42.1%, 
the  activation  loss  changes  from  30.39  V  to  41.24  V,  increasing 
about  35.7%,  the  ohmic  loss  changes  from  3.31V  to  7.83  V, 
increasing  about  136.6%. 

In  the  fuel  cell  bus  power-train,  there  is  low  voltage  protection 
function  in  the  DC/DC  converter  to  protect  the  fuel  cell  stack. 
The  protection  average  cell  voltage  is  about  0.65  V  (because 
under  the  cell  voltage  lower  than  0.65  V,  some  cells  can  not 
work  properly,  especially  in  the  large  stack  with  many  cells  in 
series.).  So  for  the  5  kW  low-pressure  PEMFC  stack,  the  output 
power  at  the  operating  voltage  52  V  was  selected  for  evaluating 
the  performance  of  the  stack.  Table  3  shows  the  performances 
of  the  stack  with  aging  under  the  protection  voltage  (52  V). 

At  the  beginning  of  the  driving  cycle  test,  the  output  power 
of  the  stack  under  the  voltage  of  52  V  is  7630  W,  and  the  current 
density  is  524.4  mA  cm-2. 

After  2500  h  driving  cycle  operation,  the  output  power  of 
the  stack  under  the  voltage  of  52  V  drops  to  2066.3  W,  and  the 
current  density  is  142.2  mA  cm-2.  The  output  power  degraded 
72.9%. 

So  under  the  driving  cycle  test  rated  current  135  A  (current 
density  is  482.14  mAcm-2),  the  stack  average  cell  voltage  will 
be  lower  than  0.65  V  after  500  h  operation.  In  the  actual  tests,  it 
did  happen  that  the  test  station  shut  down  occasionally  because 
of  one  of  the  cell  voltage  of  the  stack  lower  than  the  protected 
voltage  of  0.4  V  after  500  h  driving  cycle  test.  If  the  driving 
cycle  test  rated  current  is  90  A,  then  the  output  voltage  of  the 


Table  3 

The  performance  of  the  stack  with  age  under  the  protection  voltage  (52  V) 


Age  (h) 

Current  density 
(mA  cm-2) 

Stack  output 
power  (W) 

Power  degradation 
(%) 

BOL 

524.4 

7630 

- 

500 

417.5 

6073 

20.4 

1000 

323.7 

4706 

38.3 

1500 

247.5 

3410 

55.3 

2000 

188.3 

2735.5 

64.1 

2500 

142.2 

2066.3 

72.9 

stack  will  be  lower  than  52  V  after  1000  h  operation.  The  rated 
power  of  the  stack  must  be  reduced  if  the  stack  will  be  operated 
continuously. 

In  order  to  use  this  model  to  evaluate  the  performances  of  the 
stack  in  the  fuel  cell  bus,  the  data  of  the  fuel  cell  bus  stack  are 
compared  with  the  results  of  the  model.  The  stack  voltage  degra¬ 
dation  related  to  distance  is  almost  linear  as  shown  in  Fig.  15.  A 
set  of  stable  tests  is  carried  out  in  lab  before  the  real  route  test 
and  after  30,000  km  (about  824  h,  statistically,  the  average  speed 
of  the  bus  in  the  route  is  about  36.4  km  h-1)  test.  The  average 
cell  V-I  curves  of  the  fuel  cell  bus  stack  are  shown  in  Fig.  16. 
The  average  cell  voltage  dropped  from  0.741  V  to  0.665  V  under 
the  rated  current  90  A  after  30,000  km  test.  This  corresponds  to 
820  h  driving  cycle  tested  in  the  test  station  very  well.  It  can  be 
calculated  that  after  the  bus  operating  36,000  km  (about  989  h), 
the  average  cell  voltage  will  be  lower  than  0.65  V  under  the  rated 


Before  Route  T est  B  After  30,000  km  Route  T est 


Fig.  16.  Comparison  of  the  bus  V-I  curves  before  and  after  30,000  km  route 
test. 
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current  90  A.  This  means  that  the  output  power  of  the  stack  will 
not  meet  the  request  of  the  bus  and  the  rated  power  must  be 
reduced  if  the  bus  will  be  operated  continuously. 

6.  Summary  and  conclusions 

A  semi-experimental  voltage  degradation  model  for  a  low- 
pressure  PEMFC  stack  is  introduced  in  this  paper.  The  concen¬ 
tration  term  has  been  extracted  from  the  PEMFC  output  voltage, 
and  the  concentration  resistance  coefficient  has  been  defined. 
The  constants  of  the  model  have  been  derived  from  the  elabo¬ 
rately  designed  tests  carried  out  in  a  5  kW  low-pressure  PEMFC 
stack.  Analysis  of  the  model  shows  that  the  activation  overvolt¬ 
age  dominates  the  total  losses  with  almost  80%  followed  by 
the  ohmic  loss,  and  the  concentration  loss  is  the  mildest  and 
does  not  change  with  aging  in  the  driving  cycle  condition  (the 
concentration  loss  might  have  only  been  insignificant  due  to 
their  relatively  low  current  density  where  mass  transport  would 
not  play  a  major  role).  The  comparison  of  the  simulation  of 
the  degradation  model  with  the  actual  test  data  of  the  PEMFC 
bus  running  for  30,000  km  shows  that  the  bus  PEMFC  stack 
rated  power  or  rated  current  must  be  reduced  after  36,000  km. 
This  degradation  model  can  be  used  for  optimal  configurations, 
optimal  control  of  the  fuel  cell  system,  and  the  fuel  cell  hybrid 
power-train  system. 
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